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ABSTRACT 

Both numerical simulations and observations indicate that in an advection- 
dominated accretion flow most of the accretion material supplied at the outer 
boundary will not reach the inner boundary. Rather, they are lost via outflow. 
Previously, the influence of outflow on the dynamics of inflow is taken into account 
only by adopting a radius-dependent mass accretion rate M = Mo(r/rout)* with 
s > 0. In this paper, based on a 1.5 dimensional description to the accretion flow, 
we investigate this problem in more detail by considering the interchange of mass, 
radial and azimuthal momentum, and the energy between the outflow and inflow. 
The physical quantities of the outflow is parameterized based on our current 
understandings to the properties of outflow mainly from numerical simulations 
of accretion flows. Our results indicate that under reasonable assumptions to the 
properties of outflow, the main influence of outflow has been properly included 
by adopting M = Mo{r /routY ■ 

Subject headings: accretion, accretion disks — black hole physics — hydrody- 
namics — ISM: jets and outflows 



1. Introduction 

There are now strong observational evidences for the existence of outflow in accretion 
flow system. One of the best examples comes from Sgr A*, the supermassive black hole 
located at our Galactic center. The accretion flow in this source is likely in the form of the 
advection-dominated accretion flow (ADAF, or radiatively inefficient accretion flow; Yuan, 
Quataert & Narayan 2003). On one hand, from the observational results from Chandra 
combined with the Bondi accretion theory we can calculate the value of the mass accretion 
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rate at the outer boundary — the Bondi radius. On the other hand, radio polarization obser- 
vations constrain the accretion rate at the innermost region of the accretion flow nearly two 
orders of magnitude lower than that determined at the Bondi radius (e.g., Marrone et al. 
2006). This implies that about 99% of the material available at the Bondi radius will not 
finally enter into the black hole horizon, rather, they must be lost in the form of outfiow. 
Outfiow seems to exist also in more luminous sources whose accretion mode is different from 
the ADAF. For example, the blueshifted absorption fines, which indicates the existence of 
outflowing materials, have been detected in the X-ray spectrum of some Seyfert 1 sources 
(e.g., NGC 3783: Kaspi et al. 2001) and more spectacularly in quasars (e.g., PG 1115-1-80: 
Chartas, Brandt & Gallagher 2003). The existence of outflow has been paid more and more 
attention recently in the field of galaxy formation because of its feedback effect in the co- 
evolution of galaxy and the central active galactic nuclei (e.g.. Silk & Rees 1998; Granato et 
al. 2004; Springel et al. 2005). 

Many work has been done on the origin and dynamics of outflow (e.g., Xu & Chen 1997; 
Blandford & Begelman 2004; Xue & Wang 2005) in the frame of self-similar hydrodynamical 
solution. Magnetic field, especially its poloidal component, may presumably serve as the 
most promising mechanism on producing outflow, as proposed by, e.g., Blandford & Payne 
(1982). This has been conflrmed in non-radiative magnctohydrodynamical (MHD) numerical 
simulations of accretion flows (e.g.. Stone & Pringle 2001; Igumenshchev et al. 2003; Vlahakis 
& Konigl 2003; McKinney 2006). Radiation pressure could be another important mechanism 
in luminous accretion disk (Proga 2003). But even in the absence of magnetic held and 
strong radiation, outflow is hkely present in ADAFs. This is flrst proposed from analytical 
argument that the Bernoulli parameter of an ADAF is large or even positive because of the 
small radiative energy loss. This implies that the gas is inclined to escape once they are 
perturbed (Narayan & Yi 1994; Blandford & Begelman 1999). This suggestion was later 
confirmed by numerical simulation by Stone, Pringle & Begelman (1999). 

Since the outfiow is likely very strong (Misra & Taam 2001), they may provide an 
additional important sink of angular momentum and energy. So it is important to investigate 
its dynamical infiuence to infiow. This problem has been investigated by Kuncic & Bickncll 
(2007) in the context of the standard thin disk. In this paper we focus on ADAFs. Blandford 
& Begelman (1999, hereafter BB99) examined this question through a one-dimensional self- 
similar approach. A phenomenological way was adopted in which they parameterized the 
rate at which mass, angular momentum and energy are extracted through outflow, regardless 
the mechanism of the formation of outflow. 

BB99 gives a quite general description, covering a broad kind of outflow including Poynt- 
ing flux whose mass flux is zero while energy flux is not. It is based on self-similar assumption. 
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For the purpose of application and comparison with observation, however, we need to dis- 
cuss it based on global solutions, because self-similar solution is too simplified to be used to 
calculate the emitted spectrum. Quataert & Narayan (1999) presented the first effort on this 
aspect. They calculate the global solution of inflow when strong outflow is present by using 
a radius- dependent mass accretion rate, M oc r'^ while keeping all other equations describing 
inflow such as the momentum and energy equations unchanged. This is roughly equivalent 
to assuming that the specific angular momentum and energy of outflow is identical to the 
inflow at the same radius where outflow is launched (see e.g., eqs.([T])-([S])). This approach is 
subsequently adopted in almost all following works (e.g.. Yuan, Quataert & Narayan 2003). 
In this paper, we refer to this treatment as standard treatment. 

While the approximation of Quataert & Narayan (1999; see also Yuan, Quataert & 
Narayan 2003) may capture the most important influence of outflow to inflow, it is not 
obvious in what degree we can use this approximation or how good this approximation is 
when we compare the theoretical prediction such as the spectrum to observations. This is 
the aim of the present paper. More specifically, by considering the conservations of fluxes 
of mass, momentum, and energy of the combined inflow/outflow system, we focus on the 
influence of outflow on the dynamics of inflow. We will use a "1.5 dimension" description 
of the accretion flow, which means the height-integrated equations will be used instead of 
a fully two-dimensional description, but the conservation equations take into account the 
outflow in the vertical direction as well. The paper is organized as follows: in §2 we present 
basic equations for our model and discuss the main properties of outflow. The calculation 
results are presented in §3. The last section is devoted to a summary. 



2. Accretion Model with Outflow 

2.1. Basic Equations 

We adopt a cylinder coordinate (r, (/), z) to describe a steady axisymmetric {d/dt = 
d/d(f) = 0) accretion flow. The Paczyhski & Wiita potential (Paczyhski & Wiita 1980) 
■ip = —GMbh/{t^ — ^g) is adopted to mimic the geometry of a Schwarzschild black hole, with 
Mbh is the mass of black hole and = 2GMbh / (? is the Schwarzschild radius of the black 
hole. As shown in Fig. [1], we divide the whole accretion flow at each radius into two parts, 
i.e., inflow and outflow. For the inflow, we assume a hydrostatic balance in vertical direction 
[vz = for the inflow) and assume all quantities such as the radial and azimuthal velocity {vr 
and Vfjj), ions and electron temperature (Ti,Te) and the sound speed (c^) are only functions 
of radius r. Such an isothermal assumption in the vertical direction results in a density 
distribution of p(r, z) = p{r, 0) exp(— z^/2if^) in the inflow, where H = Cs/^Ir is the vertical 
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scale height of inflowQ. We set z = H as the surface where an outflow launches. The vertical 
gradients of above quantities are absorbed by their discontinuity between inflow and outflow 
at this surface, except that the density distribution is continuous and the density of outflow 
at z = if is then e~^/^p(r, 0). Note that the vertical velocity of outflow v^^^, will "compress" 
the inflow because of momentum conservation, thus the vertical scale height may be smaller. 
We neglect this effect here. 

From compressible Navier-Stokes equations, we can write the equations of the conserva- 
tions of mass, momentum, and energy for the inflow as follows (see Appendix A for details): 



dMjr) 
dr 



rjiAnrpv,^^ (1) 



dVr f^r,w '^r o2\ ^ dP ^ dc , , 

,d6e Pe dp, €f,yj ee„i __ , ^ 

pVr[-j ^-r)+VipVz,w 7T = dq^ + qie~q (4) 

dr dr n 

P^r{^ - + VipVz,u,^^ = (1 - S)q^ - Qie (5) 

Here all quantities have their usual meanings. The specific internal energy of electrons 
and ions are ee,ej, respectively. The pressure P is the sum of gas and magnetic pressure 
P = Pgas + Pmag- The luflow's accrctlou rate is defined as M(r) = —AnrpVrH and r]i = Pw/p 
is the density ratio of outflow and inflow (Appendix A). Parameter 6 describes the fraction 
of the turbulent energy dissipation rate q^ (= T^r^dVL/dr) that heats electrons directly. 
Energy transfers from ions to electrons through Coulomb collisions at a volume rate g^e, 
and radiative cooling rate is denoted by q~ . The quantities with subscript w denote the 
quantities of wind/outflow just away from the launching surface z = H. We take the a 
viscosity description for the stress tensor r^r (Shakura & Sunyaev 1973): 

T^r = -OlP (6) 

where a is the dimensionless viscosity parameter. Other stress tensor components are ne- 
glected for simplicity, except that the (f)z component is considered by taking into account 
the angular momentum exchange between inflow and outflow aX z = H. 



^This kind of vertical density structure is based on a first-order expand of Paczynski & Wiita potential, 
which is not exact far away from the equatorial plane (Gu & Lu 2007). 



- 5 - 



Obviously, it is impossible to directly solve the eqs. ([I])-©. We therefore introduce the 
following parameters to evaluate the radial, azimuthal, and vertical velocity, and the ion 
and electron temperatures of outflow in terms of inflow. 





(7) 


'^<j>,w = i<j)V(j,, 


(8) 




(9) 


Ti^w = ^T^Tii 


(10) 


Te,w = ^TeTe- 


(11) 



Here Vff are the free-fall velocity, v^, and Cg are azimuthal velocity and the sound speed of 
the inflow, respectively. We assume that these parameters are independent of radius. While 
this assumption is simple, we think it can capture the main physics of the influence of outflow 
in a reasonable way. Specifically, this simple assumption does not mean all quantities are a 
power-law function of radius as the usual "power-law" assumption of the mass flux of inflow. 
If we know the values of these parameters, we will be able to get the global solution of eqs. 

m-m- 



2.2. Outflow's Properties 

We now estimate the properties of outflow. Generally all these quantities should be a 
function of z. Here we consider the properties of outflow when they are just launched or 
detached from inflow. All their subsequent evolution should be due to outflow itself and does 
not affect the inflow any longer. 

The first quantity is the strength of the outflow, or the mass lost rate. BB99 assume 
M cx with < s < 1. This ensures that the mass accretion rate decreases while the 
released energy increases with accretion (BB99). The strength of outflow in our notion is 
mainly governed by C,z- The above range corresponds to < < —Vr/rjiCs ~ 0.2 (ref. eq. 
( IT2|) . but note s{r) now is a function of r). 

We next consider the value of ^<^. The vertical distribution of angular momentum of 
the accretion flow is complicated. Two-dimensional self-similar analysis on ADAF based on 
hydrodynamics shows that the specific angular momentum of outflow is lower than inflow 
(Narayan & Yi 1995; Xu & Chen 1997; Blandford & Begelman 2004). This result is con- 
firmed later by numerical simulations (e.g.. Stone, Pringle & Begelman 1999). However, any 
magnetic coupling between inflow and outflow will likely lead to transportation of angular 
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momentum from the former to latter (Spruit 1996; Stone & Pringle 2001; BB99; Bland- 
ford & Begelman 2004). Therefore in this paper we explore in a range around unity, 
0.8 < < 1.2. 

Hydrodynamical and MHD simulations also reveal that the specific internal energy or 
the temperature of the gas increases from the equator to higher altitude (e.g., Stone, Pringle 
& Begelman 1999; De Villiers et al. 2005; Beckwith, Hawley & Krolik 2008). One underlying 
reason may be that the gas with higher internal energy may escape more easily. We therefore 
consider = 1, 1.5 and = 1; 1-5. 

It is highly unclear about the radial velocity of outflow when they are just launched 
although we somehow know how they will be accelerated later. But we speculate that it 
should be positive, and should not be larger than a fraction of the local Keplerian or free-fall 
velocity Vff. Fortunately, although it may be important for the dynamics of outflow itself, 
we find that the value of has minor effect on the inflow. We simply set C,r = 0.2 in our 
calculations. 



3. Results 

For our specific model, we adopt the black hole mass Mbh = 4.0 x 10^ Mq, accretion 
rate at the outer boundary Vout = lO^r^ is M = 1.1 x 10~^MEdd, where MEdd = lOLEdd/c^, is 
the Eddington accretion rate. The values for other parameters are a = 0.1, j3 = Pgas/Ptotai = 
0.9, S = 0.3. These parameters are close to those in Yuan, Quataert & Narayan (2003) to 
model the supermassive black hole in our Galactic center. There they assume M oc with 
s = 0.27 being a constant. Under our notation, we have 

d\nM{r) Vk , . 

^^^) = ^hr;- = ^^^^^' ^'^^ 

where Vk is the Keplerian velocity. We would like to note that s{r) now is not a constant 
as in Yuan, Quataert & Narayan (2003) (or Quataert & Narayan 1999). The slope of M(r) 
now is steeper at large radius while flatter at small radius, because of the quicker increase 
of Vr compared to Vk- This is shown in Fig. [21 where we adjust the parameter C,z so that the 
accretion rates at Tout and horizon are the same as the case of M oc with s = 0.48. 

We first investigate the effect of ^z- Fig- E] shows the effects of various ^2 on inflow, with 
other outflow parameters fixed at = = ^Ti = 1-0 and Vr^^ = ^r- The four plots show the 
Mach number, profiles of density, temperature, and specific angular momentum. The dotted, 
dashed, and long-dashed lines correspond to = 0.01,0.05 and 0.15, respectively. As the 
outflow becomes stronger (^^ increases), the gas density decreases while the ion temperature 
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decreases. This is very similar to the case of the standard treatment (with increasing s). 
The decrease of ions temperature is because when more and more accretion material is 
lost via the outflow, the density profile becomes flatter thus the compression work which 
is an important heating mechanism for ions becomes weaker. Different from the ions, the 
electrons temperature has no obvious relation with the strength of outflow. This is because 
different from ions the compression work in the electron energy equation is about one order 
of magnitude smaller due to the lower electron temperature (ref. eq. Hj). 

We mentioned before that the value of the radial velocity of outflow or equivalently C,r, 
has minor effect on the dynamics of inflow. The "kick back" force due to the discrepancy of 
the radial velocity between the inflow and outflow is manifested by the second term in eq. 
([2]). Since v^^w = ^zCs, H = Cg/flk and Vr ~ av^, this term is roughly aC,z{'^ 1) times of the 
gravitational force thus can be neglected. So we simply fix C,r = 0.2 in this paper. 

We now check how good the standard treatment is. For this purpose, we first get the 
global solution with the standard treatment with M = 2 x 10~^ {r /rout)^'"^^ and Tout = lO^r^. 
We then get the global solution of eqs. ([I]) - ([S]) for various sets of outflow parameters 
of = 0.8,1.0,1.2, = 1-0 (referred to Case A) and 1.5 (referred to as Case B), and 

= 1-5. For each set of these parameters, we adjust the value of C,z so that the mass 
accretion rates at Vout and black hole horizon are equal to the values in the above standard 
treatment. By doing this, we want to focus on the influence on inflow of the transportation 
of angular momentum and internal energy between inflow and outflow, which is neglected in 
the standard treatment. Note the profile of accretion rate in this case is similar to Fig. 2. 

Fig. m shows the comparison of Case A with the standard treatment. We can see that 
our models intend to have lower densities compared to the standard treatment, although the 
accretion rates at the outer and inner boundaries are the same. We can easily understand 
this by looking at bottom-left panel of Fig. |H Our solutions have higher ion temperature 
and lower electron temperature at the inner region of the inflow. The higher ion temperature 
is because the density profile in the inner region is steeper thus the compression heating is 
stronger. The lower electron temperature is because C,t^ = 1.5 > 1 which implies that some 
internal energy is transferred into the outflow from inflow (ref. eq. S]). 

Fig. [5] shows the dynamical influences of outflow for Case B. Compared to Case A, 
the ion temperature is lower. This is obviously because C,t, is larger, C,Ti = 1-5, so some 
internal energy is transferred from inflow to outflow. The lower ion temperature results in a 
smaller H. This, combined with the smaller radial velocity, make the density of inflow higher 
compared to Case A and almost identical to that of the standard treatment, as shown in the 
figure. We also see from the figure that both the value of the specific angular momentum 
and its slope are higher compared to the standard treatment. This results in a stronger 
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viscous heating, which somehow cancel the effect of = 1-5 > 1. This is why the electron 
temperature is roughly the same with the standard treatment while higher than that of Case 
A. 

From Fig. [3] to [5l we find that within the range of the values of parameters we adopt 
to describe the outflow, the strength of the outflow has the most significant influence on the 
dynamics of inflow (i.e.. Fig. [3]). The influences of all other properties of outflow, namely 
the angular momentum, temperature, and radial velocity, are much smaller (Figs. H] and 
E]). This is the reason why the discrepancy between our model (with different properties 
of outflow) and the standard treatment, which only considers the strength of outflow but 
assume the properties of outflow are the same with inflow, is small. 

4. Summary 

Outflow is now believed to be very significant in advection-dominated accretion flow 
thus it is important to investigate their influence on the dynamics of inflow. Previously 
this was done by using a rather simple way. In the "standard treatment" the only change 
compared to the case of no outflow is that the mass accretion rate is not a constant, but a 
power-law function of radius, M oc (s > 0). All other equations describing the accretion 
flow remain unchanged (e.g., Quataert & Narayan 1999; Yuan, Quataert & Narayan 2003). 

In this paper, we investigate the influence of outflow in more detail to check how good the 
above "standard treatment" is. We have derived the height-integrated accretion equations 
including the coupling between the inflow and outflow, to investigate the influence of outflow 
on the dynamics of inflow. We assume hydrostatic equilibrium for the inflow. For the 
outflow, we assume they are launched just above the surface of the inflow. We parameterize 
and estimate the quantities of outflow in terms of the quantities of inflow mainly from the 
results of numerical simulations. In this way, we reduce the numbers of unknown quantities 
in the above inflow/outflow equations thus we are able to get their global solution. 

We have studied the influences on the dynamics of inflow of the strength (via vertical 
velocity), the (ion and electron) temperature, specific angular momentum, and the radial 
velocity of the outflow. We find that among them the strength of the outflow is the most 
important quantity. It can produce orders of magnitude difference for the density of inflow. 
If the strength of outflow is fixed, all other quantities of outflow can only produce a difference 
for the density and temperature within a factor of ~ two, if our estimations to the properties 
of outflow are correct. Therefore, the "standard treatment" is usually a good approximation. 

The largest uncertainty in our model comes from the estimations to the properties of 
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outflow such as their temperature, specific angular momentum, azimuthal and radial velocity. 
We estimate these values from numerical simulations to accretion fiows which is still not 
exact. With the rapid development of numerical simulation, our results will be significantly 
improved. Especially, if the properties of outfiow is found to be far more deviated from the 
infiow than those adopted in the present paper (e.g., >> 1), the standard treatment 
will be not enough although we believe that our conclusion that the strength of outflow is 
the most influential quantity should remain correct. In that case, the influence of the other 
properties of outflow must be taken into account as well. 

We thank our referee, Chris Done, for constructive suggestions. This work was supported 
in part by the Natural Science Foundation of China (grant 10773024), One-Hundred- Talent 
Program of China, and Shanghai Pujiang Program. 



A. The height-integrated equations describing inflow/outflow 

For a steady axisymmetric accretion flow, the equation of mass conservation is: 

Deflning the mass accretion rate M = —AirrpVrH, we integrate eq. (lAip in z from to H~^, 
here denotes just above the surface z = H where the outflow is just launched. Noting 
Vz = for inflow and Vz = Vz,w for outflow at H^, we get: 

— — = r]i4:7rrpvz,w, (A2) 

where 

P"" e~^/V(^0) n7nsn ^AQ^ 

rji = — = jj = 0.7089, (A3) 

P i/o p{r,0)exp{-^)dz 
which gives the ratio of the density of outflow and height-averaged inflow. 

The radial and azimuthal momentum equations read as follows: 

, dVr Vl dVr, dP 

P(-r^-V+-^^^ = -^ + P^r, (A4) 

PiVr'-^+^-^+Vz'-^) = \U\r). (A5) 

or r oz or 

Here Qr is the radial component of the gravitational force. The energy equation is 

pTds/dt = p {dU/dt - P/p^dp/dt) = g+ - g", (A6) 



which for steady flow reduces to: 



/dU Pdp\ fdU _Pdp\ 
\dr dr J ^ \dz dz J 



= q-q 



(A7) 



We integrate the above equations IA4tlA7l for z from z = to z = using the following 
general result: 



where / and g are functions of (r, z). Note that if g{r, z) is continuous at z = H (e.g., density 
p), the right side of equation (lASP is equal to 0. Then we will get eqs. ([2]) - ([5]). 
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Fig. 1.- 
z = H. 



Schematic diagram of inflow/outflow model 



Outflow launches from the surface 




Fig. 2. — The change of accretion rate as a function of radius for our model (sohd hne) and 
the standard treatment (dashed hne) in which M(r) oc with s being a constant. 
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Fig. 3. — Influence of the vertical velocity of the outflow described by on the Mach 
number, density, temperature, and specific angular momentum of infiow. The solid, dashed, 
and long-dashed lines are for = 0.01,0.05,0.15, respectively. The Keplerian angular 
momentum Ij. is show as the dot-dashed line in the bottom right panel. Other parameters 
are Vr,v, = Vr, = ^t, = in = 1-0. 
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Fig. 4. — Influence of tlie specific angular momentum of the outflow described by on the 
dynamics of inflow and their comparison with the standard treatment (solid line). The dot- 
ted, dashed, and long-dashed lines are for = 0.8, 1.0, 1.2, respectively. Other parameters 
are ^t- = 1 (Case A), C^j. = 0.2, C^Xe = We adjust so that all the four models have the 
same accretion rates at the inner and outer boundary. 
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Fig. 5. — Influence of tlie specific angular momentum of tfie outflow on tlie dynamics of 
inflow and tlieir comparison with tlie standard treatment. All parameters are the same as 
Fig. H (Case A) except = 1-5. 



